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Abstract. In this paper we study a linear inverse problem with a biological 
interpretation, which is modeled by a Fredholm integral equation of the first 
kind. When the kernel in the Fredholm equation is represented by step func- 
tions, we obtain identifiability, stability and reconstruction results. Further- 
more, we provide a numerical reconstruction algorithm for the kernel, whose 
main feature is that a non-regular mesh has to be used to ensure the invert- 
ibility of the matrix representing the numerical discretization of the system. 
Finally, a second identifiability result for a polynomial approximation of degree 
less than nine of the kernel is also established. 



1. Introduction 



In this work we study an integral inverse problem coming from the biology of 

►> the olfactory system. The transduction of an odor into an electrical signal is ac- 

£vq complished by a depolarising influx of ions through cyclic-nucleotide-gated (CNG) 

channels in the membrane. Those channels, that form the lateral surface of the 

cilium, are activated by adenosine 3', 5'-cyclic monophosphate (cAMP). 

• } DA. French et al. [5] propose a mathematical model for the dynamics of cAMP 

concentration, consisting of two nonlinear differential equations and a constrained 

Fredholm integral equation of first kind. The unknowns of the problem are the 

concentration of cAMP, the membrane potential and the distribution p of CNG 

• • channels along the length of a cilium. A very natural issue is whether it is possible 

to recover the distribution of CNG channels along the length of a cilium by only 

K^ measuring the electrical activity produced by the diffusion of cAMP into cilia. 

!_i A simple numerical method to obtain estimates of channels distribution is also 

proposed in [5]. Certain computations indicate that this mathematical problem is 

ill-conditioned. 

Later, DA. French & D.A. Edwards [4 studied the above inverse problem by us- 
ing perturbation techniques. A simple perturbation approximation was derived and 
used to solve the inverse problem, and to obtain estimates of the spatial distribu- 
tion of CNG ion channels. A one-dimensional computer minimization and a special 
delay iteration were used with the perturbation formulas to obtain approximate 
channel distributions in the cases of simulated and experimental data. Moreover, 
D.A. French & C.W. Groetsch [6] introduced some simplifications and approxima- 
tions in the problem, obtaining an analytical solution for the inverse problem. A 
numerical procedure was proposed for a class of integral equations suggested by 
this simplified model and numerical results were compared to laboratory data. 

l 
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In this paper we consider the linear problem proposed in [6j, with an improved 
approximation of the kernel, along with studying the identifiability, stability and 
numerical reconstruction for the corresponding inverse problem. 

The inverse problem consists in determining a function p = p(x) > from the 
measurement of 

(1) I m [p]{t) = Jo/ p{x)K m (t,x)dx, 

Jo 

for t G /, where / is a time interval, p is the channel distribution, Jo is a positive 
constant and the kernel K m (t,x) is defined by 

(2) K m (t,x) = F m (w(t,x)), 



where w(t, x), defined in (10), represents an approximation of the concentration of 
cAMP at point (t, x) and F m is a step function approximation of the Hill function 
F, given by 

(3) F(x) 



Here, the exponent n is an experimentally determined parameter and K 1/2 > is 
a constant which corresponds to the half-bulk concentration. 

Under a strong assumption about the regularity of p (namely, p is analytic), we 
obtain in Theorem 3.3 an identifiability result for (nl) with a single measurement 



of I m [p] on an arbitrary small interval around zero. The second identifiability 



result, Theorem 3.4 requires weaker regularity assumptions about p (namely, p £ 
i 2 (0, L)), but it requires the measurement of I m [p] on a large time interval. 

Furthermore, in Theorem |3.8[ using appropriate weighted norms and the Mcllin 
transform, we obtain a general stability result for the operator I m [p] for p e 
i 2 (0, L). Using a non-regular mesh for the approximation of F m , we develop a 



reconstruction procedure in Theorem 3.9 to recover p from I m . Additionally, for 



this non-regular mesh, a general stability result for a large class of norms is rigor- 
ously established in Theorem |3.11| 

On the other hand, we also investigate the same inverse problem with another 
approximation of the kernel obtained by replacing Hill's function by its Taylor 
expansion of degree m around Cq > 0. 

More precisely, the polynomial kernel approximation is defined as 

(4) PK m (t, x) = P m (c(t, x) - c ), 
where Pglji], deg (P) < m is such that 

F(x) = P m (x - co) + 0(\x - c | m+1 ), 

and c(t, x), the concentration of cAMP, is defined as the solution of the diffusion 
problem ([8]). Thus, the total current with polynomial approximation is given by 

(5) PI m [p](t)= [ p(x)PK m (t,x)dx Vt>0. 



In Theorem |8.1| we derive an identifiability result for the operator P/ m , when 
the degree of P m is less than nine. 

The paper is organized as follows. In Section 2, we set the problem, introduce 
the principal assumptions and some operator <I> m that we use to derive the main re- 
sults regarding the operator I m . These results are presented in Section 3. Section 4 
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is devoted to prove the identifiability theorems. Section 5 is devoted to the proof 
of Theorem |3.8| concerning the stability of I m . The proof of the results involving 
the reconstruction procedure are developed in Section 6, while the numerical al- 
gorithm and examples are shown in Section 7. Finally, in Section 8, we prove an 



identifiability result for PI m (Theorem 8.1 1. 



2. Setting the problem 

In this section we will set the mathematical model related to the inverse problem 
arising in olfaction experimentation. 

The starting point is the linear model introduced in [6]. As already mentioned, 
a nonlinear integral equation model was developed in [5] to determine the spatial 
distribution of ion channels along the length of frog olfactory cilia. The essential 
nonlinearity in the model arises from the binding of the channel activating ligand 
to the cyclic-nucleotide-gated ion channels as the ligand diffuses along the length 
of the cilium. We investigate a linear model for this process, in which the binding 
mechanism is neglected, leading to a particular type of linear Fredholm integral 
equation of the first kind with a diffusive kernel. The linear inverse problem consists 
in determining p = p(x) > from the measurement of 

(6) I[p](t) = Jo I P {x)K{t,x)dx 7 



dc d 2 c 

dt dx 2 

c(0,x) = 

c(t,0) = 

Yx^ L) - 


= o, 
= o, 

= c , 


t>0, x e (0,£), 

xG(0,L), 

t> 0, 


= o, 


t > 0. 



where the kernel is 

(7) K(t,x)=F(c(t,x)), 

F being given by ^ and c denoting the concentration of cAMP which is governed 
by the following diffusion boundary value problem: 



(8) 



The (unknown) function p is the ion channel density function, and c is the 
concentration of a channel activating ligand that is diffusing from left-to-right in 
a thin cylinder (the interior of the cilium) of length L with diffusivity constant 
D. I[p](t) is a given total transmembrane current, the constant Jo has units of 
current /length, and Co is the maintained concentration of cAMP at the open end 
of the cylinder (while x — L is considered as the closed end) . 

Thus, the inverse problem consists in obtaining p from the measurement of I[p] (t) 
in some time interval. 

We note that this is a Fredholm integral equation of the first kind; that is, 

(9) I\p](t)= f K(t,x)p(x)dx, 

Jo 

where K{t,x) = J F(c(t,x)) is the kernel of the operator. The associated inverse 
problem is often ill-posed. For example, if K is sufficiently smooth, the operator 
defined above is compact from L p (0,L) to L p (0,T) for 1 < p < oo. Even if the 
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operator / is injective, its inverse will not be continuous. Indeed, if / is compact 
and 7 _1 is continuous, then it follows that the identity map in L p (0, L) is compact, 
a property which is clearly false. 

In what follows, we consider a simplified version of the above problem under 
more general assumptions than those in [6] . Let us introduce the following generic 
assumptions: 

(i) We can approximate the solution c(t, x) of ^\ as follows: 

x 
(10) c(t, x) ~ w{t, x) = cq erfc ' 



where erfc is the complementary error function: 

z 

(11) erfc(z) = 1 = / cxp(-T 2 )d-r. 

V 1" J 


(ii) We consider the following approximation of Hill's function given in (pi) 

n i 

(12) F{x)^F m {x) = F{c )^a j H{x-a j ) Vxe[0,c ], 

3=1 

where H is the Heaviside unit step function, i.e. 

[ 1 if u > 0, 

(13) H(u) = \ 

[ if u < 0, 

and Oj , Oj are positive constants such that 

n i 

(14) E a i = 1 > 

3=1 

and 

(15) < a\ < OL2 < ■ ■ ■ < a m < Co, 

and hence, {aj}"L 1 defines a partition of the interval (0,co). 
With the above assumptions we define the approximate total current 

(16) I m [p](t) = Jo p(x)K m (t,x)dx, 

Jo 

where 

(17) K m (t,x) = F m (w(t,x)) = F m (c erk(-^=)) . 

Therefore, our inverse problem consists in recovering p from the measurement of 
I m \p](t) for alH> 0. 

For any 7 > 0, we consider the function <r 7 (x) = |x| 7 , and introduce the following 
weighted norms 

ll/llo,7,6 = ll Cr 7/llL 2 (0,6) ' 

ll/lli,7,6 = \W-yf\\ni(o,b) ' 

11/11-1,7,6 = 11^/11^-1(0,6)" 
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Lk = L//3k for k = 1, ...,m, and Lq = 0, 



We set 
(18) 
where 

(19) Pj = erfc _1 (a : ,/co)2V^D for j = 1, ..., m. 

On the other hand, we have 



I m [p](t) = Jo / p(x)K m (t,x)dx 



(20) 



^o^(co) /^ fflj / p(x)H(w(t, x) — atj)dx 

3 = 1 J ° 
m p 

JoF^^aj / p(x)dx, 



J ~ Gj(t)n(0,L) 

with Gj(i) = {a; G R : w(t, x) > (Xj}. Since the "erfc" function is decreasing, we 
have 



(21) 



Gj(t)= [O.&Vt], 



where {Pj}T = \ are given by (19). (Note that f3\ > /?2 > • • • > /? m .) Thus, we have 



hj(Vi) 



(22) 



/ m [p](t) = J F(c )fX;a j y 

.7 = 1 n 



p(x)dx I , 



where /ij(s) = min{L, Pjs}. 
Next, we define 



(23) 

and obtain 

(24) 
with 



$ m [<p]{t) = ^ajip{hj{t)) Vt>0, 



J m [p](t) = JoF(c )$ ro M(v/t), 



Jo 



Clearly, $ m is linear, and it follows from (14) that $ m (l) = 1, and that for any 
/ EL°°{0,L) we have 

ll*T n [/]||£« ( 0,i ra )<||/|| i oo ( 0, Z ). 

Furthermore, for any / £ C([0, L]) with /(£) = 0, we have 



(25) 



1$ 



; " ; Lno.L,,. > *-i (Z^ a i^i ' J II/IIlp(0,L) 
3=1 



1 < p < CO. 
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3. Main results 

In this section we present the main results in this paper. We begin by study- 



ing the functional $ m , defined in (23). It is worth noticing with (24 1 that the 
idcntifiability for <& m is equivalent to the identifiability for /„,.. 

Firstly, we discuss some identifiability results for the operator $ m . We begin 
with the analytic case. 

Theorem 3.1 (Identifiability for analytic functions). If tp : [0,L] — >• R is an ana- 
lytic function satisfying 

(26) * m [p](t)=0 We(0,<5) 
for some 8 > 0, then ip = in [0, L\. 

The second identifiability result requires less regularity for p, provided that a 
measurement on a sufficiently large time interval is available. 

Theorem 3.2. Let ^: [0,L]->K fen given function satisfying 

(27) *m[v»](*)=0 Vie[0,L m ]. 

Then tp = in [0, L\. 

The proof of Theorem |3.2| uses algebraic arguments and it gives us an idea on 
how the kernel could be reconstructed and also how one can envision a numerical 
algorithm. The corresponding identifiability results for the operator I m are as 
follows. 

Theorem 3.3 (Idcntifiability for analytic functions). If p : [0, L] — > K is an analytic 
function such that 

(28) I m [ P ](t)=0 V£€ (0,5), 
for some 6 > 0, then p = in [0, L\. 

Theorem 3.4. Let p : [0,L] — > R &e a ^iwen function in L 2 (Q,L) such that 

(29) / m [p](t)=0 Vt€[0,L^]. 
T/ien p = in [0, L]. 

Theorems |3.3| and |3.4| follow at once from Theorems |3.1| and |3.2| by letting 

<P(x) = / p(r)dT. 
Jo 

Let us now proceed to the continuity and stability results. 

Theorem 3.5. Let <p € H l (0,L) be a given function. Then there exists a constant 
C\ > such that 

(30) \\$m[<p]\\ m (o,L m )<Ciy\\m(o,L), 
where C\ depends only on L,j3i and fj m . 
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We are now in a position to state our first main result. Firstly, we define the 
function 



(31) 



where i 



AM = 






-(i+7-is) 



-1 is the imaginary unit. 



Theorem 3.6. Let ip £ C([0, L]) be a given function. Then there exists a constant 
7o £ R suc/i i/iai /or anj/ 7 > 7o, 



(32) 
where 



C 7 l^(-) - ^L)|| 0j%L < ||*m[v»](0 - *mM(im)|| 0j7 ,L m . 

^:=infA^( S )>0. 



It is worth noting that p2j can be viewed as an inverse inequality of ( 25 ) for 
p = 2 and for functions <p £ {/ £ C([0,L]); /(£) = 0}, and it can also be regarded 
as a stability estimate for the functional $ m . Its proof involves some properties of 
Mellin transform. Hereafter, we refer to 70 as the smallest number such that 

C 7 >0, V 7 >7 . 

Next, we present a continuity result for the operator I m . 

Theorem 3.7. Let p : [0, L] — >• R be a function in L 2 (0, L). Then, for 7 > | t/iere 
exists a positive constant C\ > swc/i i/iat 

(33) ||/ m [/=»] || lj7)iSi < Ci ||p|| ia(0ji) , 
where C\ depends only on L, cci, a m _i, a m , a m and^. 

Besides, we present a stability result for the operator I m . 

Theorem 3.8. Let p : [0, L] — > R be a function in L 2 (0,L). Then, for any 7 > 
max{7o,3/4} 7 there exists a positive constant Ci > such that 

(34) ||HLi, 7+ i,l<C*2|IUp]IIi^-I,l^ 
where C% depends only on L, C 7 > and 7. 



Theorems |3.7| and |3.8| are consequences of Theorems 3.5 and 3.6 respectively. 



Even if the proof of Theorem 3.2 is provided for any choice of the partition 
{oj'ITLi of [0, cq], its proof can be considerably simplified in the special case when 



(35) 



^ = C ° erfc I2VB ) 3= " 



,m, 



with /? € (0, 1) and (3 > constants. Note that the corresponding mesh is non- 
regular. 

In what follows, I m and $ m are denoted by I m and $ m , respectively, when oy 
is given by (35). 

For the reconstruction, we introduce the function 

/Q fi N /A ^mMft 2 //3g)-/m[p](^) w r„ „ T , 

(36) 5W = T~Eu~~\ Vt £ [0, (j L m ) . 



JoF(co) 
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As mentioned in the Introduction, we look for a reconstruction algorithm and a 
numerical scheme to recover function p from the measurement of I m [p] . We begin 
by recovering (p : [0, L] — > R, which satisfies 



(37) 



* m [0\{t/Po) = g{t), Vt<=[0,p L r , 



Next, we define functions ipi,<p2, ..., <f m by means of the following induction formu- 
lae: 



1 



(38) 



ipxix) 



0, 



, if x e l/3L,L), 



otherwise, 



and 






(39) 






ip k +i(x) = < 


a m [ g {^) g^ 1 +m {^)j 


, if ie {{3 k+1 L,{3 k L), 




°. 


otherwise, 


for k = 1 , . . , 


m — 1. Furthermore for k > m, we define 




(40) 






ip k +i(x) = < 


^^(^)-E«^^-^(p)j 


, ifaje^* 4 - 1 ^,^), 




0' 


otherwise. 



With the above definitions we have the following reconstruction result: 



Theorem 3.9. Let p be a function in C°([0,L]), let g be defined as in (36), and 
let {<Pj}j>i be given by (38)-(40). Then the function yp defined by 



+oo 



(41) 



^2ipj(x), ifx£ (0,L], 

<f>{x) = { 3=1 



5(0), ifx = 0, 

is well defined and satisfies 

(42) *.»[$(*/#)) = g{t) vt e [o, #,L m ]. 

Furthermore, p satisfies 



(43) 



P (z)dz = fts) + /m 7 [P J ( , L y Vx G [0, L] . 



Theorem 3.9 provides an explicit reconstruction procedure for both operators 
<£> m and I m and therefore a numerical algorithm for the reconstruction. 



Remark 3.10. Theorem 3.9 allows the recovery of tp, solution of (42 1, without 
any restriction about g. If another mesh is substituted to the mesh given in (35 1, 
the recovery of ip imposes to do some assumptions about g. 
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The previous reconstruction procedure gives us the possibility to obtain a sharper 
stability result. We shall provide a stability result for <& m in terms of a quite general 
norm. 

We consider a family of norms ||-||r & n for (some) functions / : [a, b) — > K, where 

< a < b < oo, that enjoys the following properties: 

(i) ||/|| M) <^forany/eiy 1 ' 1 (a,6); 
(ii) If [ai,&i) C [a, b), then 

(44) ll/ll[a lW ^IMI[a,6); 

(iii) For any A > 0, there exists a positive constant C(X) such that 

(45) ll.9A|| [Aa , Ab) <C(A) ll/H M) , 

where g\(x) = /(x/A), and C(-) is a nondecreasing function with C(l) = 1. 

A natural family of norms fulfilling (i), (ii), and (iii), is those of LP norms, where 

1 < P < +oo. Indeed, (i) and (ii) are obvious, and (iii) holds with 

C l\\ _ f A? if pe [l,+oo), 
\ 1 if p = oo. 

Another family of norms fulfilling (i), (ii), and (iii), is the family of BV- norms: 

k 

(46) ll/llw(a, 6 ) = ll/llL»(a, h )+ ^p ]T \f(x k ) - f(x k ^) \ . 

a<xi<---<Xk<b - -J 

Here, we can pick (7(A) = 1. (Note that W 1 ' 1 {a : b) C BV(a,b), see e.g. 0.) These 
kinds of norms are adapted to functions with low regularity, as e.g. step functions. 
The second main result in this paper is the following stability result. 

Theorem 3.11. Let p £ C°([0, L\) be a function and let a family of norms satis- 
fying conditions (i), (ii) and (iii). Then, we have for all k > 

(47) 

C( 8 ,n ) 
M-) - v {L)\\ [pk+1L0kL) < C(f3o)^L-L $ mM (.) _ $ m [<p](L m ) 



[/3 k + 1 L m: L m ) 



where (p(x) = / p{T)d 



Theorem 3.11 shows in particular that the value of ip in the interval [/3 fc+1 L, f3 k L) 
depends on the value of $ m [</5] in the interval [/3 k+1 L m , L m ), a property which is 
closely related to the nature of the reconstruction procedure. 

4. Proof of identifiability results 



This section is devoted to proving the identifiability results for the operator <I> I; 
Proof of Theorem 13. Ii Let ip be an analytic function such that 

n 

* m [¥>](*) = £«W(fcf(*)) - V* G (0,S). 
i=i 
Then, taking t e (0, min{5, Li}) and using the fact that 
(48) L < L-l < • • • < L m , 
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we see that hj(t) = j3jt, j — 1, ..., m. Then, we have 

m 

^ajyiPjt) = 0, te (0,mm{S,L 1 }). 

If we derive the above expression and evaluate it at zero, we obtain 

^ fe) (o) (f>^) fe j =0 Vfc>0, 

where (p( k >(0) denotes the k— th derivative of ip at zero. Since clj, j3j are positive, we 
have that YllLi a j(Pj) k > 0i therefore (p^(0) = for all k > 0, and hence (p = 0. 
This proves the identifiability for $ m in the case of analytic functions. □ 

To prove Theorem |3.2[ we need some technical lemmas. 

Lemma 4.1. Let f,g : [0, L] — > E be functions, and let s, a G [0, 1) and A € (0, 1) 
fee numbers such that 



(49) 

and 

(50) 

TTien 

(51) 

where ct\ — Amax{s,ao} 



f(r)+g(Xr) = Vte{ S L,L), 
/(r)=0 Vre[a i,L). 
flf(r) = Vre[aiI,AL), 



Lemma 4.1 is a direct consequence of (49) and (50). 

Lemma 4.2. Lei / : [0, L] — > R fee a function, and let s, ao € [0, 1) and A € (0, 1) 
&e some numbers such that 



(52) 

where 

(53) 

Then, if s > 0, 
and i/ s = 0, 



/(r)=0 Vre[5 fc L,i) Vfc > 1, 
<5fc = Amax{s, &k-i} Vfc > 1, 

/(r) = VT6[sAi,X), 



/(r)=0 Vr€(0,L). 

Proof. To prove the above lemma, we need to consider two cases: s = and s > 0. 
If s > 0, we claim that there exists fco such that ak < s. Otherwise, if a^ > 
s V/c > 0, replacing in (53), we have 

&k+i = A<5fc, 

and hence a^ = <SoA fc — > 0, which is impossible, for s > 0. 
Using (53), the desired result follows, since 

fife = As Vfc > fc . 
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Now, if s = 0, replacing it in (53) we obtain 

a k = a X ■ 



Then, using (52 1 we have 

/(t)=0, Vre(0,L), 
which completes the proof. □ 

Lemma 4.3. Let f : [0, L] — > M. be a function, and let s,a n £ [0,1), Ai,...,A„ <E 
(0, 1) and a k > 0, k — 0, ..., n be some numbers such that Ai > A2 > • • • > A„ > cxq, 
and 

n 

(54) ao/(t) + X)°i/( A i*) = ^£[sL,L), 

and 

(55) /(t)=0 VT£[a L,L). 
Then 

(56) /(r)=0 Vre[aL,L), 
where a = A„s. 

Proof. We prove this result by induction on n. 



Case ri = 1. In this case, from (54) we have the following equations 

(57) a f(t) + oi/(Ait) =0 VtS [sL, L), 

(58) /(r) = Vre[a L,L), 

and «o < Ai. Then, applying Lemma |4. 1| with g = /, we get 

/(r) = Vre[aiX,AiL), 
where ai = Ai max{s, ao}, and thus 

/(r) = VrG[aii,L), 

for a < Ai. 

If ao = 0, we obtain the desired result: 

/(r)=0 Vre[AisX,£). 



On the other hand, when ao > 0, we can apply Lemma [4.1| again with ao replaced 
by ai, since we have 

ao/(*) + Oi/(Ait) = Vt£[sL,L), 

/(r) = Vre[oiL,L), 
and ol\ < Ai. Thus, we get by induction on A: > 

(59) /(r) = Vre[a fc L,L), Vfe > 1, 
where 

(60) a k = Ai max{s, a^-i} Vfe > 1. 



12 
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Note that, if s = 0, letting t = in (|57j) yields /(0) = 0. Using Lemma |4. 2 1 with 
( 59 1-( 60 ), we conclude that 



/(r)=0 Vr6[Aisi,L), 
which completes the case n = 1. 



Case n + l. Assume the lemma proved up to the value n, and let us prove it for 
the value n + 1 . 

Assume given a function / : [0, L] — > M and some numbers s, ao G [0, 1) , a k > 
for < k < n + 1, Ai, ..., A n+ i € (0, 1) with 1 > Ai > A 2 > • ■ ■ > A rl+ i > ao, and 
such that 

n+l 

(61) a f(t) + j^ajfiXjt) = ViG[s£,Z), 

and 

(62) /(t)s0 VrG[ooi,L). 
Then we aim to prove that 

/(r)=0 VtG[A„ +1 sL,L). 
We introduce the function 

n+l , n+l 



A 



Hr) = E °i/( 1T T ) = a ^ r ) + E a i/(V)' 



i=i 



Ai 



i=2 



A, 



where Xj — ^-, j = 2, ..,n + 1. 
Ai 



Then, using (62 1, we have 
(63) 



o 



V>(r) = VrG[A 1T ^L,i). 

'Vi+l 



On the other hand, from (61 ), we have 



a /(r) + V(Air) = VtG[s£,£). 



Then, from (62 1 and Lemma 4.1 with g ~ tp, we conclude 



ip(r) =0 Vr G [Ai max{ao, s}L, AiL). 



Next, we set Si = Ai max{ao, s} G [0, 1)- Using (63), we have ^ = on [siL, Ai-L)U 
IXi^-L.L). Therefore, with ^- < i 

1 A„ + i / A„ + i — 



1+1 

n+l 



(64) 



V(r) = oi/(t) + E «i/(V) = Vr G [*ii, L) 



Note that 1 > A2 > A3 > • • • > A„+i, and that ao < A n+ i < -^ = A n +i. Then, 



by using the induction hypothesis with ( 64 ) and (1621) , we obtain 



/(r) = Vre[oiL,L), 

where a± = siA„+i = A Il+ i max{s, ao} < A n +i- Then we can repeat the latter 
argument replacing ao by &i, and we obtain 

/(r)=0 Vre[a k L,L) Vfc > 1, 
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where 

(65) a k = A n+ i max{s, a k -i} V/c > 1, 

with 0.0 = ao given. If s = 0, letting t — in (|6"T|) yields /(0) = 0. Using Lemma 



14.21 we infer that 

/(r) = 0, Vr£[aL,L), 

where a = X n+ is, which completes the proof. □ 

Proof of Theorem 13.21 Let (p : [0, L] — > M be a function such that 

771 

*m[p](*) = Yl a Mhj(t)) = Vi e [0, L m ]. 

Then, if i = L m , we obtain 

hj(L m ) = L \fj = l,...,m, 
and hence 

(66) = $ m [<p](L m ) = <p(L). 
Next, for any k £ {1, ....,m}, we have 

m 

53o^(/8j-*) = Vt€[Lfc_i,Lfc], 

which is equivalent to 

(67) a k cp(t) + J2 a W (if*) = ° W G [/3fc£fc-i,/3 fe L fe ] = [/3 fc L fc _ 1? i], 

j=fc+i \ ft / 

for k = 1, 2, ..., to. We aim to prove that 

¥>(t)=0 Vr€[/3 m L fe _i,I], 
for fc = 1, ..., to. We proceed by induction on i — to — k £ {0, ..., to — 1}. 



Case i = 0. Letting k = m in (671 yields 

am<p(t) = vt e [/? TO ir m _i, £], 

which implies 

(68) p(<r) = VT€[ft»Vi,4 

which completes the case z = 0. 



Case i = l. Letting k = m — 1 in (67), we obtain 



(69) a m -i<p(t) + a m tp ( m t ) = Vt e [/3 m _iL m _ 2 ,I]. 

,Pm-l 



We infer from Lemma 



4.3 



(applied with Ai = a Pm , s — „ 1 and an = a ) that 

v rL Pm—l' Pm-2 Pm-l' 

ip(r) =0 Vt £ [f3 m L m -2, L]. 



Case i. Assume the property satisfied for i — 1, i.e., 
(70) <p(r)=0 yr£[l3 m L m ^L}. 



14 



C. CONCA 1 , R. LECAROS 1 ' 2 , J. H. ORTEGA 1 & L. ROSIER 3 



Replacing k = m — i in ( 67 1 , we obtain 



(71) a m -np{t) + J2 a W ( -p~ l ) = ° Vi e [Pm-iL m -i-i,L] 

j—7n—i-\-l 



Pr, 



Then, if we set A 



, - - — - < 1, for j = m — i + 1, ...,m, 

Lm—i—i 



and ao 



Thus 



Pn 



S Pm—i 



A m , then we infer from Lemma 



4.3 



that 



tp( T ) = Vr e \p m L m -i-i, L\- 



tp{r)=0 Vr€[/? m L fc _i,I], 
and for fc = 1, ...,m. This implies (with fc = 1 and Lo = 0) 

<^(t) = Vr€[0,L]. 



The proof of Theorem 3.2 is complete 



□ 



5. Proofs of the stability results 



We first prove Theorem |3.5[ 
Proof of Theorem 13.51 First, some estimates are established. 



\\(p o h 



J»L 2 (0,L m ) 



1 



m ^ 2 (h 3 (t))dt = f ' V 2 {^t)dt + V 2 {L)L 

n JO \Pm Pj 

L 



(72) 



Pj JO Pwi Pm ^ J 

-^(l + ||T L || 2 i)||^||^ 1(0iL) , 



< 



where Tl(u) — u(L) is the trace operator in -ff 1 (0, L) 
Now, if we set 

1 



C\ 



lp\ 



1 + 1177, II W 



then using (72), we obtain 

(73) PmMlli^o,^) 



<^ Gj ||< y 9o/ li || L2(0Lm) <Ci||^| 



H^O.L) 



J' = l 
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On the other hand, let ip be any test function with compact support in (0,L TO ). 
Then 



$ m [<p](t)il/(t)dt = J2 



3 = 1 



, ! - ) j ip(0jt)i/j'(t)dt + (p{L) i>'{t)dt 

JLj 



(74) 



.7=1 J ° 



3=1 J ° 

where H denotes Heaviside's function. Thus 

m. 

(75) (* m fo>])'(t) = J2 a Jpjv'(l 3 M 1 - H W - £)) v * e (o,L m ). 

3=1 

Therefore, for any y> G ff 1 (0, L), the function $ m [y] belongs to i? 1 (0, £ m ). This, 
along with ( 75 1 yields 



(76) ||(*m[v])'|| ia( o lim) < 



£a.V& / (*/ 
i=i v° 



1/2 



)\t)dt) < y/&W\\„ %Il) 



Combining ( 76 ) with equation ( 73 1 , we obtain 

ll $ ™Mlli.o,i m <<5ilMli, ,£> 



where C"i = -\/(ci) 2 + /3i. The proof of Theorem 3.5 is therefore complete. □ 



Now we proceed to the proof of Theorem |3.6| Before establishing this stabil- 
ity result, we need recall well-known facts about Mellin Transform (the reader is 
referred to [12] Chapter VIII, for details). 

For any real numbers a < ft, let < a, ft > denote the open strip of complex 
numbers s = a + it (a, t £ M) such that a < a < ft. 

Definition 5.1 (Mellin transform). Let f be locally Lebesgue integrable over (0, +oo) . 
The Mellin transform of f is defined by 

M[f]{s)= f f{x)x s - 1 dx \/s^<a,ft>, 
o 
where < a, ft > is the largest open strip in which the integral converges (it is called 
the fundamental strip). 

Lemma 5.2. Let f be locally Lebesgue integrable over (0, +oo). Then the following 
properties hold true: 

(1) Let s a G R. Then for all s such that s + s Q G< a, ft >, we have 

M[f(x)}(s + s ) = M[x 8 °f(x)](s). 

(2) For any ft € R, if g(x) = f{ftx), then 

M[g}(s)=ft- s M[f}(s) Vse<a,ft>. 
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'2\ n_„ i.„„j.- .. . t2i 



Definition 5.3 (Mellin transform as operator in L ). For functions in L (0, +oo) 
we define a linear operator M. as 

M:L 2 (Q,+oo) — > L 2 (-oo,+oo), 



f 



M[f](s):=^M[f](l-is). 



Theorem 5.4 (Mellin inversion theorem). The operator M. is invertible with in 
verse 

M~ l : L 2 (-oo,+oo) — > L 2 (0,+oo), 

tp — > M" 1 ^*) := -^ f±£ x~*- iB <p(8)da. 
Furthermore, this operator is an isometry; that is, 

M[f] T2l =ll/ll^(o )0 o) V/eL 2 (0,+ M )- 

L J ( — 00,00) v y 

Proof of the Theorem 13. 6i We note that for any function / : [0, +oo[- 
that supp(/) c [0,L), we have 

f{h j {t)) = fW)- 

Thus, we obtain 



such 



(77) 



3=1 



where {/9j}yLi has been defined in (19). 

Pick any ip <G C([0, L]) and let g : [0, L m ] — > K be such that 
(78) *mb](*)=5(*) VtG[0,£ ro ]. 

Define the functions 

(79) 



?(*) = 



9 



(t)-ff(L ro ) 0<t<L. n 



y> 



•(i) - p(L) < t < L, 



t > L. 



v *>£ ro , [0 

If we replace t by L m in ( 78 ) , we have the following compatibility condition 

ip{L) =g(L m ). 
Since $ m [l] = 1, we infer that 
(80) $ m [<p](t) = g(t) Vi>0. 



Letting / = (p in ( 77 ) yields 

3=1 

It follows from Lemma fe. 21 that 

(81) X [<b m [0\] (s) = I £ ai Pj s M{0\{s) Vs E< a, (3 >, 



7 = 1 
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where < a, /3 > is the fundamental strip associated with (p. 



Let 7 > be a fixed constant. Using (81 1 and Lemma 5.2 we obtain 
(82) A£(«) M[x*v(x)](s) = M[x^ m [<p](x)](s) Vs e K, 



where A 7 ,, has been defined in (|3l|). On the other hand, 

-( 7 +i-i s ) 



A^(s) > a m /? m ' 5 

(83) > a m /3 m 7 ^ 

Therefore, if we choose 



m — 1 



E ^ 



j=i 



o-(7+i) _ o-7-t 
Hm-l ~ l Jm 



^— dmPn 



771—1 

E a A 

3=1 

(7+|)\ 



"(7+1) 



7> 



ln(a m ) 1 



then 



A£(s) > & 



(7+|) N 



Pm-1 



A, 



>0 VsG 



Thus, there exists 70 such that 

C 7 =infA^( S )>0 V 7>7 o. 



Therefore, using the fact that M. is an isometry and (82), we obtain 

(84) C7ll0llo, 7f L<ll*4flllo l7 ,^. 

which completes the proof of Theorem |3.6| 



D 

We are now in a position to prove Theorems |3.7| and |3.8| 

Proof of Theorem 13. 7i Let us fix any 7 > and let p : [0, L] — > K be a function 
in L 2 (0, L). From p4| we have 

(J^WW)' = 7^ 7 " 1 /77 l M(x) + ^(/,„[p](x))' 
(85) = 7x 7 ' 1 /7»M(x) + X ^" J ° F{C0) (*mM)'(V5), 

X 

where ip(x) — J p(r)dr. (Note that (p € -ff 1 (0, L).) Since 


x^ 1 {^ m [p}Y(Vx)) 2 dx = 2 /^ r 4 ^ 1 (($ m M)'(r)) 2 dT = 2 W^WWl^^ 

JoF(c )\ 



we have 

ll^mWIIl.V,^ < ll I m[p]|lo, 7 ,Z,= l + (7ll J m[p]|| , 7 -l,L^+ " 



l^mWWo^t,! 



< H^HIo,^ + 2 7 2 l|/ m Wllo, 7 -i,^ + ( J ° F ( C °)) 2 IK $ -M)'llo,2 7 -i 



(86) < (i 2 + 27 2 )||/7„M|| 2 i7 _ 1 , L?n + (Jo^( C o)) 2 ||(<i> m M)'llo, 27 -i 
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On other hand, using (24 1 and the change of variable t — x 2 , we have 

II^MIllU-l,^ = , F , \ T ,2 f^ ^-" (lm[p]( X 2 )) 2 dx 

2 (F{co)Jo) Jo 

1 



2(F(c )J 



2 ll 7 m[P]llo, 7 -l,I4 



By replacing (87) in (861, we obtain 
\\I m [p}\\l^ L 2 m < (L 2 +2 7 2 )2(^( Co )Jo) 2 ||$ m M||^ 27 _3 !Lm + (F( Co )Jo) 2 ||(<i> m M)'||o !27 _i iLm , 

we have 



and assuming that 7 > |, from Theorem 
III, 



3.5 



m^jII 1>7i zs, < V^ +4 7 2J F( C0 )L 2 ^I ||$ m M|| H1(0 , Lm) 

< v/3i 2 + 4 7 2 Jo^(co)L 2 ^-iC 1 || V || H1(0>£) . 
But, from Cauchy-Schwarz inequality we have |<^(x)| < \L ||p|| L2 / L \ , and hence 

o n o n o 

IMIffi(o,£) = IMIl=(o,l) + ||p||z,2 (0 ,ri <(L +1) ||p|| L 2 (0 ^) . 
Therefore, for any 7 > f , we have 



\l 



m\P\ Nl,7,££, 



< Ci Hp|li2(o,z,) 



where 



d = V3£ 2 +4 7 2 Joi 71 (c )£ 27 - i C , i(£ 2 + 1) 1/2 , 
and the proof of Theorem |3.7| is therefore finished. 



□ 



Proof of Theorem 13.81 Let ip be any test function compactly supported in (0, L), 
and let 7 be a positive constant. Set 



g 1 {x) = x 7 p(x), 



tp(x) = / p{r)dr 
Jo 



and 

It follows that 

and hence, 

< g-y+i,i> > 



<p(t) = <p{x) - <p(L). 

(x 1+1 fi(x))' = (7 + l)x J 0(x) + g 1+ i{x), 

L i-L 

g 1+ i{x)ip(x)dx = I ((a; 7+1 (^(x))' — (7 + l)x 1 (p(x)) tp{x)dx 



(x~ f+1 Lp(x)ij'(x) + (7 + 1)^ 'ip(x)ip(x)) dx. 
jo 

Then, we have 

|<5 7+ i,V>| < (||^IIo, 7 +i,l + (7+1)II^IIo, 7 ,l) MhHqx) 



Therefore, 

(89) 



< (£ + 7 + l)Hv|lo, 7 ,LllV'l| ff i(o,i)- 
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Thus, using Theorem 3.6 we have that for any 7 > max{7o, |} there exists a 
constant CL > such that 



IMI-1,7+1,£ - ll37+lllff-i(0,L) 



(90) 



T I 



<(L + 7 + l)C- 1 ||$ ro MI| 0>7jim + ^jf=\ 



$ m [v?](L m )| 



Using ( 24 ) , we have 
(91) 



<f> m [<p](L m ) 



F(c )J { 



■Im\p](L 2 m ). 



Replacing (91) in ( |90| ) and using (871, with 27 — 3/2 replaced by 7, we obtain 
(1 + 7 + 1) 



||p||_i, 7+ i,i < 
Therefore, setting 



(7-1 J 1 Jo ' 
72|JoF(c )| 7 I V27TT 1 



II^IIMI^WHi. 



7 _ 1 r 2 



(£ + 7+1) ^-1 f -, /q I'm 



\T, 



L'i, 



a 



we obtain (34). The proof of Theorem 3.8 is achieved. 

6. Numerical reconstruction results 
This section is devoted to the proof of Theorems |3.9| and |3.11| 

Proof of Theorem 13.91 Let p be a function in C°([0, L]), and let us consider the 
functions {ipj}j>i defined in (38)-(40). 

First, we note that for all k > 1 we have 

(92) ¥>*(&) = 0, Vx^^L^^L). 

Then, we can define the sequence {ipp} p£ iq<- as 

p 

tpp{x) = 2, <Pj( x ) ^ x £ R- 

Using ([92]) we have that for all X £ R \ (0, L), 

ij) p (x) = VpeN*, 

and hence, 

lim ijiJx) = VxeR\ (0,L). 

p— >+oo 

Besides that, we consider the ceiling function 

\x] = minjfc G Z \ k > x}, 
i.e., \x~\ is the smallest integer not less than x. 
Next, we define 

~ln(x/L)~ 



(93) 

Then, we have 



k*(x) 



Vse (0,L). 



ln(/3) 
x e [^'^L,^'^- 1 ^, Vx € (0,L). 
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Therefore, we obtain for x € (0, L) 

ip p (x) = <Pk*(x)(%) Vp>k*(x), 
and hence, 



(94) 



lim tj) p (x) = ip k * {x) {x) \/x e (0,L). 



Thus, the series in (41 1 is convergent, i.e. the function (p is well defined. 



On the other hand, by replacing (351 in (191 we obtain 



By using ( 92 ) we have 



ft =/3o0 i , i = l,..-,m. 



<p(x) = Q, V.Tel\[0,I), 



and combining with ( 77 ) , we get 
(95) 



*m[#(t/A)) = $>iW*) Vt > 0. 



By replacing i = and £ = /3 L TO in (95), and using (|36|, (38), and (41 ), we obtain 



<t> m [<p](0) = g(O), and $ m [$(L m ) = g(/3 L m ) = 0. 
Now, if we take t £ (0,/3 £m) = (0,L//3 m ), we have 

Fte(p,0*- m L), for je{l,2,...,m}. 
We need to consider two cases, t < L and t > L. 
Case t < L. In this case we have 

/^G(0,i), for je{l,2,.,m}, 
and 



Thus, replacing (94) in (95), we obtain 



$m[$(*/A>) =53a i ^ i+fc .(t)(/8 J '*), 



and hence, using (40) with k + 1 = m + k*(i), we obtain 



$m [$(*/$)) = J! ^fj+fc'W^fj + ttm^+t^tjtr') 



m— 1 



771—1 



J=l \ .7=1 

= .9(0- 



Case t > L. Let us set 



*♦(«) 



ln(x/L) 

>(!//?). 



Va;>L, 
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where 

[x\ = max{fc 6Z k < x}, 

is the floor function, i.e., it is the largest integer not greater than x. Thus, we have 

(96) /3 k *( x)+l x<L</3 k * {x) X Vx>L, 
and 

(97) k*(/3 k * {t)+ H) = l. 



Then, we infer from ( 96 1 that 

0H>L Vj<fc*(t), 

0H<L Vj>h(t) + l. 



By using (95) and (94), it follows that 

771 

*m[0](*/A>) = E a Wk*{PH) (fl J t) 

j=fe,(t)+l 
m—k* (t) 
= E a 3+k,{t) l fk*(fii + k -Wt) [/3 : ' +fc * (t) iJ 



From (97), we have 



k*(pi+ k »V>t)=j Vj>l, 



and hence, from ( 38 )-( 39 ), with k + 1 = m — fe*(t), we obtain 

m—k* (t) 

* m [m/fa) = E *j+*.(t)^ (>' +Mt) *) 



m— fc«(t)-l 

E a i+ ; - 

m-fc«(t)-l 



.(t)Vj (>' +fe * (t) *) +o m v m -*.(t) (/3 m *) 



i=i 



-fe.(*)-i 

E 



.9(0 - E «M*)+J^ f>' +Mt) *) 



<?(*)■ 



It remains to prove ( 43 ) . Replacing (|36j) in ( |42| and using ( 24 ) , we get 

/ m [p](i 2 //3o 2 )-£n[p]04) 



*m[fl(*/A>) 



J F(c ) 



(98) = $ m M(V/3o)- / "; [ y^ ) vte[o,A,x m ], 



where ip(x) = J p(r)dr. Using $ m [l] = 1 and Theorem 3.2 we obtain (43), i.e. 

i m [p]{Ll " 



(p{x) = <p(x) ~ 



JoF(co) 



V.t e [0, L] 
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This completes the proof of Theorem |3.9[ 



□ 



Proof of Theorem 13. Ill Let p be a function in C°([0, L]), and let {<fij}j>i be 
defined in @-(|40|. Using @, we obtain 

(99) ip(x) = ip(x) - ifi(L) \/x e [0, L], 

where tp — L p(r)d,T and (p has been defined in (41 ). 

Recall that the family of norms || • ||[ a ,b) (for < a < b < oo) satisfies (44)- (45 1. 



Using (|41|, (|93]-(|94| and ((99j), we obtain 

(!00) \\<p(') - V>(L)\\ [0k +i L ^k L) = ll^fc+ill^+iL^L) 

Let us prove that for any k > 0, we have 



(101) 



C{(3 m ) „ 

Wk+lW^k + iL.pkL) < k+1 \\9\\[p^ + ip L m ,l3 L m ) 



The proof of ( 101 ) is done by induction on k. 



Case k = 0. Using (38) and (44)-(45), we have 

C((3 m ) „ 

WVlWlfi L,L) > — \\9\\[pi3 L m ,P L m ) > 

as desired. 

Assume now that for all j = 1, ..., k (with k > 1), we have 



(102) 



1 1 1 1 \P J 1 1 II 

WiWlfliLJJ-i-L) - ~j \\9\\[pip L m ,poL m ) ' 



and let us prove (101). 



Case k + 1 < m. Using (39) and (44)- (45 1, we obtain 

ll^/c+lll^fc+iL.^L) - ~~~ (C(/3 m ) \\9\\[pk+i p L m ,p k p L m ) 



v- r(P k+1 \\i II 



i=i 



Using induction hypothesis (102), we have 

IIVM-l||[£*+i£,/9 fc £) ^ 



-i- C(/3™)||. 9 || [ , fc+1/3oLm ^, oLm) +^a„ l _ fc _ 1+J c(^) ^|| fl || 

- 1 ™ \ 7 = 1 \ P J Ctm 



[pi P L m ,p Q L m ) 



< 



C(/3 



- Q k+T I a m \\9\\[p^ip L m ,P"p L m ) + 2-^ a m-k-l+jC f -^- J ||#|| \pi p a L m> p oLm ) 









i=i 
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Note that 

(103) C{u)<\ Vue(0, 1), 

for C(-) is nondecreasing and C(l) = 1. Therefore, 



C 



P j 



<i Vie {i,..., fc}. 



Thus, we obtain 



ii 1 1 \M / 1 1 II 

Wk+lW^pk+i^k^ < fe+ i ll3ll[^fc + i/3 L m ,/3 L m ) 



This proves (101 1 for all k = {0, ...,m - 1}. 



Case fc + 1 > m. Replacing </?fc+i by its expression in (40) and using (44|-(45) 
and the induction hypothesis, we obtain 



1 / m_1 / B m \ 



B m \ C(0 m ) „ 



fe+1 °m llfl'll[ ; efc+l / 8 i m ,^* / 9 D i m ) + E °^ 

Am \ =1 



^-y nyii^+^^+i/So^^oL™) 



- „k+l "m 



J=l 



/T 



2_ ^ ^^J 1 IIS , ll[^ S! + 1 /5oim,/3oi m ) ' 



the proof of (101 1. 



and with (103 1 we infer that C(/3 m /^) < 1 for all j £ {1, ..., to- 1}. This completes 



On the other hand, using ( 42 ) and (W5| , we obtain 

ll5ll[/3fc + i / g L m ,/3 L m ) ^ C (M 






[/3 fc + 1 L m ,L m ) 



By replacing (99) in (101), we obtain 



C(B m ) 
ll¥>k+l||p*+iL,/9*£) ^ W-y iM(') 



[^ + ii m) i m ) 



and by replacing in ( 100 ) , we obtain ( 47 ) . This completes the proof of Theorem 3.11 



□ 



7. Numerical results 



In this section we discuss the numerical implementation of the scheme developed 
when proving Theorem |3.9| 
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Firstly, we define {a}" l =1 as in (35), and let 



Fi = < 



p i ocj + a J+l 



j = l,...,m- 1, 



, F ( c o) j = m, 

where F is Hill's function defined in pi). Next, we set 

(104) aj = Fj ~ Fj -\ j = l,...,m. 

r(c ) 

Since F is increasing and < F(x) < 1 for all x > 0, we infer that a,j > for all 
j = 1, ..., m, and that 

m 
3=1 

The corresponding approximation F m of Hill's function is shown graphically in 
Figure [T] 




Figure 1. Hill's function and its approximation 



Recall that a non-regular mesh in the interval [0, L] was introduced when proving 
Theorem |3.9[ Now, let us start defining 

Pq,i = {(xo,Xi,...,x q ) E M. q+1 : Xj E [0L,L), x = f3L, xj-i < Xj, Vj = 1, ...,q} , 

and its representative vector 

P 1 = (xo,x l7 ...,x q )ER c * +1 . 
Next, introduce the sets 

V q ,, = {x\ p j - l x£V q ,x), ]>l 
and denote their corresponding representative vectors by 

Pj = P 3 ~ 1 V l = /3Pj_i e R q+1 . 
Let us fix some p > 1. Our aim is to recover the function p on the mesh 

(105) £ p , g =uf =1 p<w, 
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where the corresponding representative vector is given by 

P = (P l! P 2 ,...,P p )er+'+ 1 . 



By using (38)-(40), we can define the vectors Gi, G2, ..., G m £ R q+1 inductively as 
follows: 



(106) 

and for k = 1 , . . . , m — 1 



(Gl)s = ^Hl^ 



, 8 = 1,. 



1, 



(107) (G fc +i) s = ^ ffl(^ i )-E°-*-i+i(Gi).j 
Finally, we define the vectors G& € M 9+1 for k = m, *..,p — 1, by 

m — 1 



•5=1, 



,g + i. 



(108) (G 



fc+ljs — 



Introduce the vector 



9 



k+l)s 



0T, 



~ Z^ °i (Gj+fc-m+l) s 
3=1 



!,...,«+!. 



G = (G 1 ,G 2 ,- 



•,G p )e 



? p+g+i 



which represents a discretization of the function (p given by Theorem |3.9| on the 
mesh defined by £ Pl9 ,; that is, ((P) s , (G) s )^_' is a discretization of the curve 
(x,<p(x)), x E (0,L). 

Therefore, using ( 43 1 and applying a forward difference scheme, we obtain an 
approximation of the curve (x, p{x)), x e (0, L), through the vectors X, Y € R p+q 
given by 

(109) (X) S = (P) S , (Yj^m^l gj^-g^ O 



s = l, 



,P- 



It should be noted that the maximum function was considered in ( 109 ) because of 
the positivity restriction on the density function. 



7.1. Examples. Let us consider 

8a 8 x 7 
(110) p{x) = , 

(x s + apy 

with a = 1.5. Figure u\ shows functions p{x) and ip(x) defined in (110) and their 
approximations obtained by the previous procedure. 



<P(x) = / p{r)dT = 
Jo 



x 

x s + a» : 



Remark 7.1. If we consider any discretization of (16) on a given mesh, one has 
to solve a system like 

Ay = g. 

Obviously, the system depends strongly on the choice of the mesh. 

We notice that the matrix A may not be invertible. While it is difficult to give a 
general criterion for the invertibility of A in terms of the mesh, Theorem \3 . S\ guar- 
antees that the matrix A is indeed invertible when the non-regular mesh described 
in (105) is used for the discretization of system (16). 
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A: Function p. 
B: Function <p. 



" 


J\ U;r... 


■ 




- 


/ \ 




( f i 1 




Figure 2. The target functions p and ip with their approximations. 
Let us now consider the example studied in [BJ. To this end, we dehne 

0, t G (0,t£) e i a y), 



(111) 



I(t) 



1 



£ ^Delay 



t > i-Delay 



with tpeigy — 30 [ms], n/ 



2.2, J Maa! = 150[pA] and Kj ~ 100[ms]. The current 
given in ( 111 ) is a sigmoidal function with short delay (Figure wp), which is similar 



to the profiles encountered in some practical situations (see e.g., [I], [3] or pTj). 

The numerical solution corresponding to these data is shown in Figure [3]A.. It 
should be noted that the numerical solution given here is perfectly consistent with 
those obtained in [3 . 

8. Polynomial approximation of Hill's function 



In this section we consider the same inverse problem with another approximation 
of the kernel in fcf\ , for which we keep the function c and replace Hill's function F 
by a polynomial approximation around Cq. 

More precisely, let P m be the standard Taylor polynomial expansion of degree 
ttt, of ([3J around Co; that is, P m € M.[X], deg(P m ) < m and 

(112) F(x) = P m (x - c ) + 0(\x - c \ m+1 ). 
A new approximation for the kernel is defined by 

(113) PK m (t,x)=P m (c(t,x)-co), 
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A: approximation of p. 



B: current I(t) as defined by (fill 







1 \ 




J 1 2 3 4 5 




200 400 600 300 1000 1200 1400 



Figure 3. Approximation of the function p(x) with current I(t) 
as defined by (fff ) 



where c(t, x) is the solution of Q, given by 
(114) 



2 \ i /+2? e -"* D ' 
c(t,x) = co - c ( T 1 2^ — ip k (x) ) , 



with 
(115) 
and 

(116) 



Mfc 



ipk(x) 



\k=0 

2k + 1 

7T. 

1/2 



/'A- 



sin(/x fe x). 



Besides, we define the total current associated with this polynomial approximation 
as follows: 



(117) 



PI m [p](t) = p(x)PK m (t,x)dx Vt > 0. 
Jo 
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Next, we present our main result regarding the operator PIm] it asserts that the 
identifiability for the operator PI m holds when m < 8. 

Theorem 8.1. Let m < 8 be a given integer. Then 

Ker PI m = {0}, 

where 

Ker PI m = {/ e L 2 (0, L)| PI m [f](t) = Vt > o} . 



8.1. Proof of Theorem 8.1, Let us start noting that {ipk}k>o is an orthonormal 
basis in L 2 (0, L). Thus, for any / € £ 2 (0, L), we can write 



(118) 
where 



f(x) =J2<f^k> M*) in L 2 (0,L), 



k>o 



</,V>fc >- 



f(x)ipk(x)dx. 



We write 



P m (z) = a + aiz + • • • + a m z m , 
where a.j £ R, for all j = 0, 1, ..., m, and introduce the set 



(119) 






Let e > be a given positive constant. For any p £ L 2 (0, L) and for all s > 0, 
we have 

PI m [p](e + s) = a / p(a;)da; 



1 = 1 "^o \fc>0 



,-nlD{e + s) 



3 = 1 

do I p(x)dx 
Jo 



E a ^~ Co 

3 = 1 



/'■A- 



-^ fe (a;) p(x)dx 



E 



E^l^( s + £ ) 



p=i 



fe 1 ,...,fe,->0 



n^^M),^,. 



o P=1 Mfc 



Note that the convergence of the last series is fully justified, as for any j e {l,...,m} 
and any s > 0, we have 



E 



-E^l^ s+£ ) 



P =i 



fci,...,&»>0 



Jo 



j V'fcp 0*0 



A*fc p 



)p(x)dx 



< 



\P\\L l (0,L) 



E 

, fe>0 



2 e -^ De 
Pk 



< oo 



Therefore, there is a family {aA(e,p)}AeA 7Tl such that 
(120) ^ |a A (£,p)|<oo 



AeA„ 
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and 



(121) PI m [p](s + s) = a / p(x)dx + V a x {e,p)e XDs Vs > 

- 70 AGA m 

Lemma 8.2. Let p e L (0, L) 6e a given junction, such that 

(122) PI m [p](t)=0 Vi>0. 
TTien 



(123) 

and 

(124) 



p{x)dx = 



a A (e,p)=0 VAeA., 



Proof. First, the series in (121) is uniformly convergent for s > 0, by (120) 
Define {AA;}fc>i as 

(125) Ai =min|A 6 A m |, A fc +i = minJA 6 A m \{Ai,..., A*}}, 

and note that this defines an increasing sequence < Ai < A2 < • • • 
Using definition 125 we can rewrite (|121[) as 



(126) PI m [p](e + s) = a I p{x)dx + y Va Xk {e,p)e- XhDe e- XkDs Vs > 0. 

Jo fe>i 



Set 



Then we have that 



-A fc De„ — AfcDs 



Vs>0. 



(127) 



|S,-(*)|<e- A ^'[X;K(e,p)|e- A * 



A fc D £ 



Plugging (|126j into ([122]), it follows that 

(128) ao/ /?(«)<& + Si (a) = Vs>0. 



Passing to the limit as s — > +00 in ( 



Si(s) -> by (1271, we obtain (123 



128), and noting that a$ — F(c ) ^ and that 



The proof of ( 124 ) is done by induction on j > 1. 



Case j=l. Plugging (123) in (128) and multiplying by e XlDs , we obtain 
(129) a Xl (e,p)e- XlDe + e XlDs S 2 {s) = Vs > 0. 



But, using (127), we also have that 

|e AlDs S' 2 (s)| < Ce-l x '- Xl)D '. 



Thus, noting that Ai < A2 and passing to the limit as s — > 00 in (129), we obtain 



.-ill 
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Case j=n+l. Assume that 
(130) a A ,(e,p) = Vj = {!,..., n}. 



Plugging |130| ) and ([123]) in (J122J) , we infer that 
(131) 



a Xn+1 (e, p)e~^+^ s + e^+^ s S n+2 (s) = 0, V.s > 0. 



On the other hand, using ( |127 ), we have that 

\e^ lDs S n+2 {s)\ < c e -(A„ +2 -A n+1 )D S _ 



Thus, noting that A„+i < X n +2 and passing to the limit as s — > oo in (131), we 
infer that 



a A„ +1 (e,p) = 0. 



This yields (124). This complete the proof of Lemma 8.2 



□ 



Lemma 8.3. Let {/ife}fc>o be the sequence defined in (115). Assume that 

.2 



(132) 



Mni H + /4 fc = M« 



for some k > 1 and n,n\, ...,rik > 0. XTiera 
(133) fc = 1 mod 8. 

Proof. We have 

7T 2 

where </?(n) = n 2 + n. Thus, substituting this expression of // 2 . in (132) yields 

k 



Noticing that <p(n) is an even number for all ngff, we obtain (133) 



D 



wo 



Proof of the Theorem 18.11 Let p £ L 2 (0, L) be a given function such that ( 122 ) 
holds. From Lemma [8. 2| we infer that ( |124| holds. If m < 8, using Lemma 8.3 
have that for all n > 0, all k £ {2, ..., m}, and all m, ...,Tifc > 0, 

2 , ,2/2 



Then, with A = /j 2 G A m , we obtain 



<V 2 (£,P) = e M " D£ ai(-co\/T) 



i>n(x) 



L Jo Mn 

Since ai = F'(co) ^ (F being increasing), we infer that 

< ip n ,p>=0, Vn> 0, 



p(x)dx. 



and hence, with (1181, that 



P = 0. 
This completes the proof of Theorem |8.1| 
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